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Abstract 

Using 3.1 fb _1 of data accumulated at the T(4S) by the CLEO-II detector, corresponding 
to 3.3 million BB pairs, we have searched for the color-suppressed B hadronic decay processes: 
B° —>■ D°(D*°)X°, where X° is a light neutral meson 7r°, p°, rj, r\ or cu. The D*° mesons are 
reconstructed in D*° — ► D°7r° and the D° mesons in D° — > K~ty + , K~ty + ty and K~n + n + iT~ 
decay modes. No obvious signal is observed. We set 90% C.L. upper limits on these modes, 
varying from 1.2 xl0~ 4 for B° -> D°7r° to 1.9 xlO" 3 for B° -> D*°r]' . 

PACS numbers: 13.25.Hw, 13.25.-k, 13.30.Eg, 14.40.Nd 
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I. INTRODUCTION 



The B hadronic decays B° — > D (D*°)X , where X° is a light neutral meson 7r°, p , 77, rj or 
have not yet been observed. These decays proceed via the internal spectator diagram shown in Figure 
^ The internal spectator decays are expected to be suppressed relative to the decays that proceed via 

2030897-001 




Figure 1: Internal spectator diagram of B hadronic decays B° — ► D°(D*°)X°. 

external spectator diagrams, since the color of the quarks from the virtual W must match the color 
of the c quark and the accompanying spectator antiquark. Therefore these decays are refered to as 
color-suppressed decays, while decays via external spectator diagrams are refered to as color-favored 
decays. Measurements of these color-suppressed decays allow tests of the factorization hypothesis 
and provide useful information on the scale of strong final-state interaction in the B meson system. 

Previous CLEO papers [2| reported upper limits on these color-suppressed B hadronic decays. 
Here we present new results using the full CLEO-II data set and an improved analysis method. 

II. DATA SAMPLE and EVENT SELECTION 

The data used in this analysis were produced in e + e~ annihilations at the Cornell Electron Storage 
Ring (CESR) and collected with the CLEO-II detector [3]. The intergrated luminosity is 3.1 fb™ 1 at 
the T(4S) resonance, which corresponds to (3.32±0.07) x 10 6 BB pairs, and 1.6 fb -1 at energies just 
below BB threshold (henceforth referred to as the continuum). 

Hadronic events are selected by requiring at least three charged tracks, a total detected energy 
of at least 0.15 E c . m , and a primary vertex within 5.0 cm along the beam {z) axis of the interaction 
point. To suppress continuum background, we require that the ratio of second to zeroth Fox- Wolfram 
moments R2 [U determined using charged tracks and unmatched neutral showers be less than 0.3 (0.5 
for clean decay modes involving i] or rj ). To further reduce continuum background, we then require 
that the cosine of the angle between the sphericity axis of the B meson candidate and the sphericity 
axis of the remainder of the event satisfies \cos(9 sp h er icity)\ < 0.8 (0.9 for decay modes involving r\ or 
i] ). For a jet-like continuum event, the two axes are almost parallel, while they are almost uncorrelated 
for a BB event, 
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III. B RECONSTRUCTION 



Selection of D° and D*° candidates 

The D° candidates are reconstructed in the decay modes D° — > K~Tt + , K~ir + ir° and K~ir + 7r + -K~ 
(charge-conjugate modes are implied). The 7r° candidates are formed by combining two showers whose 
invariant mass is within 2.5a of the 7r° mass (where henceforth a denotes r.m.s. resolution). Charged 
tracks are required to be consistent with coming from the interaction region in both the r — <f> and r — z 
planes. The measured specific inoization (dE/dx) of charged kaon and pion candidates are required 
to be consistent to within 2a for kaon candidates and 3<r for pion candidates. Charged tracks are 
required to have a momentum greater than 250 MeV for D° — > K~ir + and D° — > K~tt + tt° candidates 
and 200 MeV for D° — > K~tt + it + -k~ candidates. For D° — > K~ir + ir° decay mode, we select regions of 
the Dalitz plot with large amplitude to further suppress the combinatoric backgrounds. The invariant 
mass of D° candidates is required to be within 2. Oct of the nominal D° mass. 

The D*° candidates are reconstructed using decay mode D*° — > D°7r°. We form D*° candidates 
by D° candidates using above selection, then require that the D*° — D° mass difference be within 
2.5a of the measured value 

Selection of the light neutral meson X° 

We reconstruct 7r° candidates as described previously. The p° candidates are reconstructed in the 
mode p° — > 7r + 7r~. 

Candidate r] and r/ mesons are reconstructed in their r\ — > 77 and 77 — > r]ir + ir~ decay modes. The 
absolute value of the 77 decay angle is required to be less than 0.85 to remove asymmetric candidates 
which are primarily background. The invariant mass of each rj and r/ candidate must be within 30 
MeV of their nominal mass. 

The lo mesons are reconstructed in the decay mode u> — > ir + ir~ir°. Charged and neutral pions are 
required to have momenta greater than 250 MeV, to reject soft pions from D*° or D* + decays. The 
00 candidates are also required to be within 30 MeV of the nominal lo mass. 

All charged pion candidates used in X° reconstruction are required to have a measured dE/dx 
within 3a of the expected value for pions. 

Selection of the B candidates 

The L> ( * )0 candidates are combined with a light X° to form a B meson. At CLEO the energy of 
the B meson is the same as the beam energy and the measured beam energy is more precise than the 
reconstructed B meson energy. Full reconstruction of B mesons at CLEO makes use of this fact by 

defining two variables. One is the beam-constrained mass, Mb = \J E% eam — P observed- ^ ne other one 
is the difference between the reconstructed energy and the beam energy, AE = E b servec i — Eb eam . The 
AE variable is sensitive to missing or extra particles in the B decay, as well as particle species. For 
fully-reconstructed B meson decays, the Mg distribution peaks at 5.28 GeV with resolution around 
2.7 MeV, and AE peaks at 0.0 GeV with a resolution ranging from 18 to 50 MeV, depending on the 
B and D° decay modes. 

Since signal and background are in general much better seperated in AE than in Mb, instead of 
cutting on the AE variable and fitting Mg as in previous analyses, we cut on Mg and fit the AE 
distribution for the signal yield. 
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IV. BACKGROUND STUDY 



In our search for the color-suppressed B hadronic decay modes B° — > D^ (tt , p°, r], rj , uj), there 
are backgrounds to these decays from continuum and BB events. The continuum backgrounds are 
suppressed using event-shape variables. They are not expected to show any structure in the AE dis- 
tributions. The 1.6 fb _1 continuum data set is used to monitor the continuum background levels. We 
find the continuum background level to be very low for all color-suppressed modes. No accumulation 
around AE = is observed in the continuum data. 

The backgrounds from BB events are dominated by feedthrough from color-favored two-body 
hadronic decays of the type: 

B~ ^ D°(7T-,p-,a^), B~ ^ D*°(ir-,p-,ai) 
B° D+(ir-,p-,ai), B° -» L»*+(7r-,p-,ar) 

The branching ratios of these color-favored B meson decay modes were measured previously by CLEO 
[2]. In most cases the background arises when a real, energetic D° or D*° from the 2-body color-favored 
decays is combined with a fake light meson. 

The backgrounds from these color- favored processes can have structure in the Mb and AE dis- 
tributions, depending on which color-suppressed mode is being analyzed. Particularly important are 
color-favored B meson decays that give exactly the same final state particles as our color-suppressed 
signals do. Neither misidentification nor additional particles are needed for those color-favored decays 
to fake some signal modes. Therefore, the Mb distribution from these physics background peaks at 
5.28 GeV while its AE distribution peaks at 0.0 GeV, exactly as the color-suppressed signal. While 
D°ir° is not susceptible to this background D°p° and D°uj are, as shown below. 

Color - suppressed : B° D°p° -> D°tt + tt- B° -» D q uj D q tt + tt-tt° 
Color - favored : B° -> D*+ir- -> D°it + it- B° D* + p~ D°-k + tt--k° 

Another background that can show structure is color-favored decay in which one of the final state 
particles is lost. Examples include: 

Color - suppressed : B° L> vr° B° -» D°p° D°tt + tt- 

Color - favored : B~ D°p~ L> vr° (tt - ) B° -» D* + p~ -» £)°7r+7r- (tt°) 

These background events can peak in around 5.28 GeV when the missing 7r~ or ir° from the p~ 
decay is very soft and does not contribute much to the beam-constrained mass calculation. However, 
the AE for these background events differs from zero by more than one pion mass, due to the missing 
7r~ or 7T° from the p~ decay. For these types of color-favored backgrounds, the color-suppressed signals 
are much better separated from background in AE. 

For decay modes involving rj or r] , combinatoric background is the dominating source. Therefore, 
backgrounds for these color-suppressed processes have no accumulation in the Mb and AE distribu- 
tions. 

For B° — > D*°X°, there is no corresponding color- favored B meson decay that fakes our signal as 
B° — > D* + tt~ fakes B° — ► D°p°. Also the background level from color- favored B meson decays is very 
low for B° —s- Z?*°X° decay processes, due to the good resolution on the D*° — D° mass difference. 

Almost all the discrimination power against color-favored physics backgrounds come from selection 
cuts on X°. We make full use of mass, momentum, decay angle and other kinematic variables of X° 
to suppress backgrounds while keeping signal efficiency as high as possible. 

The X° candidates in B° -> L>M°X are very energetic due to the hard spectrum of two-body B 
meson decays. We require the momentum of the ir° candidate to range from 2.1 GeV to 2.5 GeV. 
Similar momentum requirements are imposed on the other light neutral meson X° candidates. 
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For B — > D p decays, there are color-favored physics backgrounds from B — > Z)* + 7r~ that give 
exactly the same final state particles. The B~ — > decay can also fake our color-suppressed signal 

by substituting the soft tt° from p~ decay by a soft tt + from the other £> meson. In these physics 
backgrounds, the ir~ is always much more energetic than the n + from D* + — > D°7r + decay. There 
exists a correlation between the D° and the fast 7r~ (slow 7r + ) from the fake p°. To suppress these 
physics backgrounds, we require that the D° candidate to be associated with a fast tt + (slow 7r~) from 
the p° candidate. There is still a contribution from color-favored physics backgrounds even after this 
requirement, because a D° decay has a certain chance of being misidentified as a D° decay. For the 
D°'s from our signal process, together with the dE/dx and D° mass requirements, this misidentification 
rate is determined to be less than 20%. After further suppression due to the p° mass and momentum 
requirements, the contribution from color-favored physics background is negligible. Since the p° from 
B° — > D°p° decay is longitudinally polarized, we also cut on the p° decay angle (the angle between 
the direction of the pion in the p° rest frame and the direction of the p° in the lab frame) to reduce 
combinatoric backgrounds. 

Signal selection efficiencies for all the color-suppressed decay modes are shown in Table 1. The 
systematic error due to the detection of charged and neutral tracks, together with the Monte Carlo 
statistical error, are included in the error on the efficiency for each decay mode. 

Table 1: Selection efficiencies and yields of all color-suppressed modes. The three efficiencies 
and yields of each B° -> D°(D*°)X° correspond to the three D° -> K~ir + , R-n+n and 
K~7t + 7t + 7t~ modes. 



Decay Mode 


Selection Efficiency 


Yield 


B {> -> D°7T U 


26.1±2.2, 7.8±1.0, I2.5±f.3% 


-0.3±6.4, -6.7±4.3, -3.3±7.0 




14.1±1.8, 3.7±0.7, 5.4±0.9% 


2.5±2.6, 5.0 ±3.4, -1.2±3.4 


B° -> D°p° 


8.4±0.4, 2.6±0.3, 3.9±0.3% 


1.4±3.0, -3.0±4.3, 3.1±5.0 




4.0±0.4, 1.0±0.2, I.5±0.2% 


-1.0±1.4, 1.4±1 .6, 0.8±1.3 


B° -> D°rj 


24.5±3.0, 7.0±1.2, 11.4±1.6% 


-1.4±2.0, -3.4 ±3.1, -6.0±4.0 


B° -> D*°r] 


10.5±1.8, 3.4±0.8, 4.9±0.9% 


0, 0, 


B° -> D°r)' 


13.4±1.9, 3.6±0.7, 5.9±1.0% 


0, 0.8±2.2, 1.8±3.0 


B° -> D*°rj' 


5.9±1.1, 1.7±0.4, 2.5±0.5% 


0, 0, 1 


B {) -> D°c^ 


12.4±1.3, 2.8,±0.4, 3.2±0.5% 


-4.1±4.0, 6.2±3.8, 3.6±5.6 




4.8±0.7, 1.0±0.2, 1.3±0.2% 


1.8±1.2, 0.8 ±1.8, -0.2±1.2 



V. RESULTS 

The Ai? distributions for the on T(4S) and continuum data samples of all the color-suppressed 
signal processes after all cuts are shown in Fig. 2 — 6. The AE distribution of each color-suppressed 
mode is fit with a Gaussian and a background shape function. The mean value and width of the 
Gaussian distribution are fixed with values determined from signal Monte Carlo. We use various 
color-favored decay modes B~ -» D°tt-, B~ -» D°p~, B~ -» D* -r- , B -» D*°p-, B° D+tt", 
B° — > D + p~ to check that the AE resolutions in data and Monte Carlo agree well. Possible differences 
between data and Monte Carlo in the AE distributions are considered and included in the yield error 
as systematic errors. 
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Various AE background shape functions have been used to fit for the signal yield: a simple second- 
order polynomial; or a background shape using Monte Carlo simulation BB events plus a continuum 
component represented by a second-order polynomial. For the latter shape, the BB contribution is 
scaled to the known luminosity while the continuum component is allowed to float. Our results are 
found to be insensitive to different background shapes, and both of them describe the AE distributions 
reasonably well. Differences in the yield due to the choice of AE background shape are included in 
the yield error to account for the systematic uncertainties. For each signal process with several D° 
decay submodes, the yield for each D° submode is obtained seperately, since the AE resolutions are 
different for the different modes. The results are shown in Table 1. The yields of the D° submodes 
are added independently to get the total yield. 

The formulas used to calculate the branching fractions are: 

Br (B° -» D°X°) = , ^ = (1) 

N BS x (Ef =1 Efficiency(i) x Br(D°i)) x l\Br(X°) 



Br(B° -> D*°X°) = (2) 

N BB x Br(D*° -> D°ir°) x (£f =1 Efficiency^) x Br(D°i)) x Y\Br(X°) 

where N G 6 S is the total yield summed over the three D° submodes, N BB is the number of BB pairs, 
Efficiency(i) is the selection efficiency for B° — > D°(D*°)X° decay in the ith D° submode, Br(D°i) 
is the branching ratio of the ith D° decay mode, and Y\Br(X°) is the product over all the relevant 
branching fractions of the X° decay chain. Particle Data Group values for D°, D*°, rj, rj and uo 
branching ratios are used in the upper limits calculation [5j and are listed in Table 2. 



Table 2: Particle Data Group branching ratios that are used in the upper limit calculation for 
color-suppressed B hadronic decays. 



Decay Mode 


PDG Branching Ratio 


D° -> K--K+ 


4.01±0.14% 


D° -> K-n + n° 


13.8±1.0% 




8.1±0.5% 




63.6±2.8% 


P° — > 7T + 7T~ 


100% 


rj -y 77 


38.8±0.5% 




43.7±1.5% 


uo — > 7r + 7r _ 7r° 


88.8±0.7% 



The upper limits of color-suppressed branching ratios are determined by the method described 
in section 17 of the Particle Data Group 90% C.L. upper limits on branching ratios of color- 
suppressed B hadronic decay processes, together with theoretical predictions j^], are shown in Table 
3. Among all the decay modes, the upper limit for the B° — > D°tt° mode is the lowest at 1.2 xl0~ 4 . All 
the upper limits on branching ratios are still higher than theoretical predictions Ej ■ Compared with 
factorization and QCD based calculations, no dramatic enhancement of color-suppressed B hadronic 
decay branching ratios is observed, indicating that there is no sign of large scale final-state interaction 
in these B meson decay modes. 
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Table 3: 90% CL. upper limits in branching ratios of all color-suppressed modes, together with 
comparison with theoretical predictions. 



Decay Mode 


Branching Ratio (@90% CL.) 


Theoretical Predictions 


B° -> D°7r u 


< 1.2xl0~* 


0.7xl0~ 4 




< 4.4xl0~ 4 


l.OxlO" 4 


B° -> DV 


< 3.9xl0~ 4 


0.7xl0" 4 




< 5.6xl0" 4 


1.7x10-* 




< 1.3x10-* 


0.5xl0" 4 


S° -> D* u 7] 


< 2.6xl0~ 4 


0.6xl0" 4 


B° -> I? ?/' 


< 9.4xl0~ 4 




B° -> 


< 19xl0~ 4 




£° -> 


< 5.1 xlO" 4 


0.7xl0" 4 


B° -> D* u ^ 


< 7.4xl0" 4 


1.7x10-* 
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Figure 2: AE distributions of B° — > D n° and B° — > D*°7r° decay modes. Solid histograms 
are the AE 1 distributions of the 3.1 fb _1 of data collected on the T(4S) resonance, which are 
fit using background and signal functions. Dashed histograms are from the 1.6 fb _1 continuum 
data sample. 
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Figure 3: AE distributions of B° — > D°p° and 5° — > _D*°p° decay modes. Solid histograms 
are the AE distributions of the 3.1 fb _1 of data collected on the T(4S) resonance, which are 
fit using background and signal functions. Dashed histograms are from the 1.6 fb _1 continuum 
data sample. 
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Figure 4: AE distributions of B° — > D ?? and B° — > D*°7] decay modes. Solid histograms are 
the A.E distributions of the 3.1 fb _1 of data collected on the T(4S) resonance, which are fit 
using background and signal functions. Dashed histograms are from the 1.6 fb~ x continuum 
data sample. 
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Figure 5: AE distributions of B° — > .D // and 5° — > decay modes. Solid histograms are 

the distributions of the 3.1 fb _1 of data collected on the T(4S) resonance, which are fit 
using background and signal functions. Dashed histograms are from the 1.6 fb _1 continuum 
data sample. 
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Figure 6: AE distributions of B° — * D°u and 5° — > D* u decay modes. Solid histograms are 
the AE distributions of the 3.1 fb _1 of data collected on the T(4S) resonance, which are fit 
using background and signal functions. Dashed histograms are from the 1.6 fb _1 continuum 
data sample. 
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